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ABSTRACT 


The  attitude  of  a satellite  refers  to  the  rotational  orien- 
tation of  the  spacecraft  relative  to  some  reference  triad  of 
Cartesian  axes  (these  being,  for  the  type  of  spacecraft  treated 
here,  the  orbit  radius  vector,  the  normal-to-the-orblt  plane,  and 
the  vector  cross  product  of  the  two).  Mathematically,  the  atti- 
tude Is  usually  represented  by  nine  direction  cosines  and/or 
three  Euler  angles.  The  numerical  determination  of  these  parame- 
ters Is  the  objective  of  attitude  estimation.  Various  schemes 
have  been  developed  and  used  by  the  Applied  Physics  Laboratory  to 
determine  the  attitude  performance  of  Its  satellites.  In  recent 
years,  a least-squares  technique  that  Involves  eigenvalue  and 
eigenvector  computation  has  been  added.  This  report  presents  the 
formulation  of  the  technique  and  discusses  Its  successful  applica- 
tion. Attitude  estimation  results  from  three  orbiting  spacecraft 
are  Included. 
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1.  INTRODUCTION 


This  report  presents  the  formulation  and  application  of  a 
least-squares  attitude  estimation  technique  for  gravity-gradlent- 
stablllzed  spacecraft.  The  following  topics  are  included: 

(a)  definition  of  satellite  attitude,  (b)  statement  of  the  atti- 
tude estimation  problem,  (c)  least-squares  solution  (see  Appen- 
dix A for  derivation) , (d)  computational  algorithm  for  solution 
implementation  (including  the  auxiliary  equations  derived  in  Ap- 
pendix B) , (e)  sources  of  estimation  error,  and  (f)  results  from 
several  orbiting  APL  satellites. 
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2.  ATTITUDE  ESTIMATION  PROBLEM 

GRAVITY-GRADIENT  STABILIZATION 


A gravity-gradient-stabilized  satellite  has  one  of  Its  axes 
(usually  the  Z-axls)  always  pointed  toward  the  earth.  Such  a 
spacecraft  Is  designed  to  take  advantage  of  the  fact  that  the 
earth's  gravitational  field  will  tend  to  stabilize  a trlaxlal  body 
(one  with  unequal  principal  moments  of  Inertia),  with  Its  principal 
axis  of  minimum  Inertia  aligned  with  the  local  vertical  and  Its 
axis  of  maximum  Inertia  aligned  with  the  normal-to-the-orbit  plane 
(see  Fig.  1)  (Refs.  1,  2,  and  3).  (The  local  vertical  is  an  Imagi- 
nary line  from  the  earth's  mass  center  to  the  satellite's  mass 
center. ) 

Spacecraft  built  by  APL  have  used  extendible  booms  to 
achieve  a favorable  moment -of -Inert la  distribution,  i.e.  , an  In- 
ertia ellipsoid  where  the  smallest  principal  Inertia  is  at  least 
an  order  of  magnitude  less  than  the  others.  The  satellites  dis- 
cussed In  this  report  have  also  Included  a constant-speed  rotor 
with  Its  spin  axis  aligned  with  (or,  in  some  cases,  defining)  the 
spacecraft  Y-axls.  The  addition  of  the  wheel  enhances  the  overall 
stabilization  by  adding  gyroscopic  stiffness  and  stability  to  the 
alignment  of  the  Y-axls  (Ref.  4) . 


SPACECRAFT  ATTITUDE 


The  attitude  of  a satellite  refers  to  the  rotational  orien- 
tation of  the  satellite  axes  relative  to  some  reference  triad  of 
Cartesian  axes.  For  a gravlty-gradient-stablllzed  spacecraft,  this 

Ref.  R.  E.  Flschell,  "Magnetic  and  Gravity  Attitude  Sta- 

blllzatlon  of  Earth  Satellites,"  ARS  J. , Vol.  31,  September  1961. 

Ref.  2.  R.  A.  Nidley,  "Gravitational  Torque  on  a Satellite 
of  Arbitrary  Shape,"  ARS  J. , Vol.  30,  No.  2,  1960. 

Ref.  3.  R.  E.  Roberson,  "Gravitational  Torque  on  a Satel- 
lite Vehicle,"  J.  Franklin  Inst.,  Vol.  265,  January  1958. 

Ref.  4.  V.  L.  Plsacane,  "Three-Axis  Stabilization  of  a 
Dumbbell  Satellite  by  a Small  Constant-Speed  Rotor,"  APL/JHU 
TG  855,  October  1966. 
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Fig.  1 Stabilization  of  a Satellite  with  Three  Unequal  Principal 
Moments  of  Inertia 
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reference  system  Is  called  the  local  vertical  system  (Z^  is  the 
outbound  local  vertical,  Is  the  normal-to-the-orblt  plane,  and 
Xj  Is  the  vector  that  completes  the  right-hand  set) . 


The  orientation  (or  attitude)  of  one  reference  system  to 
another  can  be  mathematically  represented  by  a direction  cosine 
(or  orthogonal  transformation)  matrix.  This  nine-element  represen- 
tation can  be  reduced  to  three  by  expressing  the  transformation 
matrix  as  a product  of  three  single-axis  rotations.  The  argument 
of  each  rotation  is  called  a Euler  rotation  angle.  Although  the 
Euler  angles  are  not  unique,  they  do  provide  a more  meaningful  and 
convenient  representation  of  the  attitude  of  one  reference  system 
relative  to  another  (Refs.  5 and  6). 


Three  Euler  rotation  angles,  R,  P,  and  Y (corresponding  to 
roll,  pitch,  and  yaw  angles,  respectively),  are  used  in  this  re- 
port to  define  a satellite's  attitude  relative  to  the  local  verti- 
cal frame.  Figure  2 shows  the  part  each  angle  plays  in  the  total 
definition  of  the  spacecraft's  orientation.  The  sequence  of 
single-axis  angular  rotations  from  which  Fig.  2 is  obtained  is  as 
follows : 


1.  The  first  rotation,  R2(P)f  is  about  Y^; 

2.  A rotation,  Ri(R),  about  the  new  X-axls  is  performed, 
and 


3.  A rotation,  R3(Y) , about  the  new  Z-axls  is  performed. 
The  axes  resulting  from  the  final  rotation  are  the 
fixed  vehicles  axes,  3^,  Yy,  and  Z^- 


The  product  of  the  three  rotations  is  the  direction  cosine  matrix, 

A,  l.e. , 


A - g^(Y)  R^(R)  R2(P) 
or 


(1) 


/R  P Y + 

P Y 

R Y 

R P Y 

- P Y 

18  8 

c c 

c s 

s c s 

S • 

A -I  R P Y - 

P Y 

R Y 

R P Y 

+ P Y 

“ 1 s s c 

C 8 

c c 

see 

8 1 

\r  P 

-R 

R P 

' C 8 

8 

c c 

Ref.  5.  H.  Goldstein,  Classical  Mechanics, 
Publishing  Co. , Inc. , Reading,  MA,  1950. 


Addlson-Wesley 


Ref.  6.  G.  A.  Smith,  "Four  Methods  of  Attitude  Determina- 
tion for  Spin-Stabilized  Spacecraft  with  Applications  and  Compara- 
tive Results,"  NASA,  TR  R-445,  August  1975. 
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where  the  subscripts  s and  c denote  the  trigonometric  sine  and 
cosine  functions,  respectively.  By  definition,  A is  the  orthogonal 
transformation  matrix  from  the  local  vertical  reference  system  to 

the  satellite  reference  system.  I 

The  functional  relationship  between  A and  the  attitude 
angles  is  very  important.  It  provides  the  connection  between  the  j 

least-squares  solution  and  the  attitude  angles.  Thus,  given  the 
computed  matrix.  A,  the  three  attitude  angles  can  be  computed  as 
follows: 

I 

Roll  - sin  ^(-A  ,) 

-1  ' 

Pitch  ■ tan  (^31^^33^ 

Yaw  - tan  ' 


PROBLEM  STATEMENT 


I 


The  primary  objective  of  attitude  estimation  Is  the  deter- 
mination of  A.  This  unique  representation  for  the  attitude  can 
then  be  used  (as  shown  above)  to  obtain  the  more  convenient  and 
meaningful,  but  nonunique,  Euler  rotation  angles. 

In  general,  the  satellite  attitude  estimation  problem 
usually  poses  the  following  question:  Given  a vehicle  with  fixed 
reference  axes  rotated  relative  to  a set  of  '»  -own  reference  axes, 
what  is  the  direction  cosine  (or  orthogonal/  aatrlx  that  relates 
the  two  reference  systems?  Or,  given  a set  of  vectors  (m*,  mj, 

...  , m**)  relative  to  the  fixed  vehicle  axis  and  a corresponding 

“Tl 

set  (m, , m_,  ...  , m ) relative  to  the  inertial  axes,  what  Is  the 
— 1 —L  “n 

orthogonal  matrix.  A,  that  satisfies  the  following  equation 


m ) - A(m  , m,,  ...  , « ) ? 


(3) 


In  most  satellite  application's,  the  first  vector  set  is  usually  the 
output  resulting  from  measurements  by  on-board  sensors  (star 
trackers,  sun  sensors,  etc.),  and  the  second  vector  set  is  the 
known  counterpart  of  the  first  set. 
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3.  LEAST-SQUARES  SOLUTION 


MATHEMATICAL  EQUATIONS 


The  least-squares  approach  to  the  problem  can  be  stated  as 

follows:  Given  the  set  of  vectors  (m*  ml,,  ...  , m*)  and  (m,  , m., , 

— 1 — n —1  —2 

...  , m^)  as  defined  previously  (for  n i 2) , find  the  orthogonal 

matrix.  A,  that  brings  the  first  set  Into  the  best  least-squares 
coincidence  with  the  second  set.  That  Is,  find  an  A that  mini- 
mizes the  scalar 


Q(A) 


E 


T 


j-1 


A T ^ 

(nij  “ 


(4) 


where  e.  Is  the  column  vector  of  errors  associated  with  the  m*  ob- 
-J  -n 

served  vector  and  the  superscript  T denotes  transposition. 

This  least-squares  problem  has  been  treated  In  the  open 
literature  (Refs.  7,  8,  and  9).  In  Appendix  A the  derivation 
(taken  from  Ref.  7)  of  a closed  form  solution  for  A Is  discussed. 
The  solution  is 


.¥.  “1  q.  ^ 

(P  ) (P  P) 


(5) 


Ref.  7.  J.  L.  Farrell  and  J.  C.  Stuelpnagel,  "A  Least 
Squares  Estimate  of  Satellite  Attitude,"  Problem  6501,  SIAM  Rev.. 
Vol.  8,  No.  3,  July  1965. 

Ref.  8.  P.  B.  Davenport,  "A  Vector  Approach  to  the  Algebra 
of  Rotations  with  Applications,"  NASA,  TN  D-4696,  1968. 

Ref.  9.  L.  Fralture,  "A  Least-Squares  Estimate  of  the 
Attitude  of  a Satellite,"  AIAA  J.  Spacecr.  Rockets.  Vol.  7,  No.  5, 
May  1970. 
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where  the  matrix,  P,  Is  defined  as 


* T 

P - M M 


(6) 


fV 

and  the  matrices,  M and  M,  consist  of  the  juxtaposed  column  vec 
tors  m*,  ...  , m*  and  m_ , respectively.  For  the  case 


and  m^ , 


respectively.  For  the  case 


In  which  the  determinant  of  P Is  less  than  zero,  the  solution  for 
A Is 


A 


(P^)  (P^P)^  (I  - 2g'^H  G) 

'■I  ' m m mm  m' 


(7) 


where  G Is  the  model  matrix  whose  columns  are  the  eigenvectors  of 
the  matrix  (P^P) , and  H is  an  nth-order  matrix  with  all  elements 
equal  to  zero  except  the  (n,  n)  element,  which  is  equal  to  one. 
The  matrix,  I,  is  the  nth-order  identity  matrix. 

The  matrix,  (P^P)^,  is  the  square  root  of  the  symmetric 
matrix,  (P^P) , with  positive  eigenvalues.  A computational  algo- 
rithm for  determining  this  matrix  contains  the  following  steps: 


1.  Find  the  eigenvalues  and  the  normalized  eigenvectors 
of  the  matrix  (P^P) ; 

2.  Construct  a matrix,  G,  in  which  the  kth  column  is  the 

* * 

kth  eigenvector  associated  with  the  kth  eigenvalue; 

3.  Construct  a diagonal  matrix,  D,  in  which  the  (k,  k) 
element  is  the  square  root  of  the  absolute  value  of 
the  kth  eigenvalue;  and 

T 

4.  Compute  the  desired  matrix,  (P  P)  , according  to  the 
equation 


T T 

(P  P)  - G D G . 


(8) 
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4.  COMPUTATIONAL  ALGORITHM 


The  actual  detennlnation  of  A Involves  a series  of  computa- 
tional steps.  For  the  applications  discussed  here,  the  following 
steps  were  required  to  Implement  the  least-squares  technique: 

* * * 

1.  Compute  the  first  set  of  vectors,  (m,  , m.,,  ...  , m ), 

— i ~n 

In  satellite  coordinates.  The  first  set  of  vectors  Is  provided  by 
attitude  sensors  on  board  the  spacecraft.  Two  types  of  sensors,  a 
triad  of  orthogonal  vector  magnetometers  (Fig.  3)  and  digital  solar 
attitude  detectors  (DSAD's)  (Fig.  4),  were  used  on  the  gravlty- 
gradlent-stablllzed  satellites  discussed  later.  The  geomagnetic 
field  vector  (denoted  m*)  and  the  sunllne  vector  (denoted  m*)  are 

obtained  from  the  output  signals  of  the  magnetometers  and  sun  sen- 
sors, respectively.  A third  Independent  vector,  m*,  can  be  com- 
puted as  the  vector  cross  product  of  the  other  two7 

2.  Compute  the  corresponding  set  of  vectors,  (m^^,  ^2*  ...  , 

m ) In  local  vertical  coordinates.  The  vector,  m, , Is  obtained 

from  the  evaluation  of  a complex  mathematical  model  (48-term  spheri- 
cal harmonic  expansion)  of  the  earth's  magnetic  field  (Ref.  10). 

The  computation  of  the  sunllne  vector,  m2.  Is  based  upon  the  cata- 
loged ephemerls  of  the  sun  (Ref.  11).  As  In  step  1,  the  cross  prod- 
uct of  and  ^ Is  used  to  provide  m^.  The  satellite's  orbit  Is  a 

required  Input  In  these  computations  because  the  mathematical  formu- 
lations are  referenced  to  an  Inertial  coordinate  system  called  the 
geocentric  reference  system  (Z  Is  the  North  Pole,  X Is  the  first 
point  of  Arles,  and  Y Is  the  vector  cross  of  Z and  X).  In  Appen- 
dix B,  the  transformation  of  vectors  from  the  geocentric  system  to 
the  local  vertical  system  Is  discussed. 

o * 

3.  Construct  the  matrices,  M and  M.  Two  3 by  3 matrices 
are  constructed  using  the  two  computed  vector  sets.  The  columns 
of  the  matrix,  M*,  are  composed  of  the  Juxtaposed  vectors,  m^, 


Raf.  10.  J.  C.  Cain  et  al.,  "Computation  of  Che  Main  Geo- 
magnetlc  Field  from  Spherical  Harmonic  Expansions,"  NASA/GSFC, 
NSSDC  68-11.  Greenbelt,  MD,  May  1968. 

Ref.  11.  The  American  Ephemerls  and  Nautical  Almanac,  U.S. 
Government  Printing  Office,  Washington,  DC,  1972-1977. 
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and  m*.  The  matrix,  M,  la  similarly  constructed,  using  the  vec- 
tors , ^ m^ , and 

4.  Compute  the  orthogonal  transformation  matrix,  A.  The 
actual  computation  of  ^ Is  preceded  first  by  the  computation  of 
the  matrix,  P (Eq.  6)  and  then  by  the  required  eigenvalue  and 
eigenvector  computations  (see  Ref.  12  for  a computational  algo- 
rithm). The  condition,  determinant  of  P 0,  that  must  be  satis- 
fied Is  violated  whenever  the  geomagnetic  field  vector  and  the 
sunllne  vector  are  parallel  (or  antiparallel).  However,  this 
rarely  occurs  when  these  two  vectors  are  used  (Ref.  13). 

5.  Compute  roll,  pitch,  and  yaw  rotation  angles.  The 
appropriate  elements  of  A are  used,  according  to  Eq.  2,  to  com- 
pute the  desired  Euler  rotation  angles. 

A con9>uter  program  was  written  to  execute  the  above  steps. 
The  satellite  attitude  results  that  are  discussed  later  were  out- 
puts from  that  program. 


Ref.  12.  S.  S.  Kuo,  Computer  Applications  of  Numerical 
Methods . Addlson-Wesley  Publishing  Co.,  Inc.,  Reading,  MA,  1972. 

Ref.  13.  H.  D.  Black  et  al. , "Attitude  Determination  Utiliz- 
ing Redundant  Sensors,"  Proc.  4th  Intern.  Aerospace  Instrumentation 
Symp. , Cranfleld,  England,  March  1966. 
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5.  SOURCES  OF  ERROR 


The  types  of  errors  that  can  degrade  the  accuracy  of  the 
least-squares  estimate  include  errors  associated  with  the  attitude 
sensors.  Inherent  in  the  least-squares  theory  is  the  assumption 
that  the  errors  are  Gaussian  (Ref.  14),  l.e.,  random  noise  with 
zero  mean.  This  means  that  the  technique  is  generally  most  accu- 
rate for  errors  that  can  be  characterized  as  Gaussian  in  their 
statistics. 

However,  systematic  errors  (e.g.,  sensor  biases),  seldom 
qualify  as  Gaussian  in  their  statistics.  Thus  this  type  of  error 
poses  the  greater  threat  to  the  accuracy  of  a least-squares  esti- 
mate. 


Although  there  are  many  sources  of  such  errors  that  can 
theoretically  degrade  the  accuracy  of  the  estimates  of  the  satel- 
lite's attitude,  some  of  these  errors  are  expected  to  be  signifi- 
cant, while  others  are  not.  Thus,  a thorough  discussion  of  all 
the  errors  will  not  serve  a useful  purpose  in  this  report.  In- 
stead, a brief  discussion  of  the  anticipated  significant  error 
sources  is  presented. 


MAGNETOMETER  ALIGNMENT  UNCERTAINTY 


There  are  always  small  errors,  called  uncertainties,  accom- 
panying the  positioning  or  aligning  of  a sensor.  An  uncertainty 
in  this  application  is  the  angle  between  the  actual  and  assumed 
position  vectors  of  a sensor.  Theoretically,  the  three  magnetome- 
ters will  be  aligned  with  the  X,  Y,  and  Z satellite  reference  axes 
so  that  the  magnetic  field  components  along  these  axes  can  be  mea- 
sured. However,  in  practice,  each  sensor  can  be  positioned  only 
to  within  0.6“  of  the  desired  vehicle  axis  (Ref.  15). 


Ref.  14.  R.  Deutsch,  Estimation  Theory,  Prentice-Hall,  Inc., 
IIJ,  1965. 

Ref.  15.  Private  communication  with  B.  Tossman,  AFL,  March 
1971. 
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SUN  SENSOR  ALIGNMENT  UNCERTAINTY 


Each  of  the  DSAD's  can  be  mounted  to  within  0.1*  of  its 
assumed  position  vector  (Ref.  16).  The  results  of  the  error  are 
analogous,  in  description,  to  those  attributed  to  the  magnetometer 
alignment  uncertainty.  However,  the  magnitude  of  the  DSAD  align- 
ment uncertainty  is  so  small  that  its  effect  on  the  attitude  esti- 
mates will  be  insignificant. 


MAGNETOMETER  BIAS 


Biases  on  the  vector  magnetometers  could  result  from  using 
incorrect  calibration  tables  or  from  actual  residual  magnetic 
dipoles  on  the  spacecraft.  Although  prelaunch  procedures  are  de- 
signed to  minimize  the  problem,  postlaunch  schemes  such  as  recall- 
bratlon  or  bias  estimation  are  usually  available. 


MAGNETOMETER  NOISE 


Because  of  sensor  electronics  nonlinearity  or  other  Inherent 
limitations,  an  uncertainty  called  noise  accompanies  any  sensor 
measurement.  For  the  vector  magnetometers,  the  expected  noise 
level  is  about  1 mA/m  (Ref.  17).  This  is  usually  described  as 
Gaussian  noise,  with  zero  mean  and  a standard  deviation  of  1 mA/m. 

The  effects  of  these  errors  on  attitude  estimation  results 
can  be  examined  via  digital  computer  simulations.  Table  1 sum- 
marizes the  results  of  a study  of  the  effects  of  magnetometer 
noise  and  alignment  uncertainty.  Three  simulations.  Cases  I,  II, 
and  III,  were  performed.  Each  covered  a time  period  greater  than 
one  orbit  (^  100  min).  Least-squares  estimates  of  the  attitude 
angles  were  computed  every  2 min. 

Case  I dealt  with  the  effects  of  magnetometer  alignment  un- 
certainty only.  The  results  of  the  investigation  indicated  that  a 


Ref.  16.  Private  communication  with  G.  Fountain,  APL,  March 
1971. 

Ref.  17.  Private  communication  with  F.  Mobley,  APL,  March 
1971. 
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Table  1 

Results  of  Testing  Least-Squares  Technique  with 
Corrupted  Magnetometer  Data 


Case  I 

Case  II 

Case  III 

Statistics 

Roll 

Pitch 

Yaw 

Roll 

Pitch 

Yaw 

Roll 

Pitch 

Yaw 

Average  Absolute 
Deviation 
(deg) 

0.07 

0.18 

0.07 

0.21 

0.30 

0.15 

0.41 

0.78 

0.39 

Standard  Deviation 
(deg) 

0.04 

0.04 

0.03 

0.14 

0.19 

0.13 

0.33 

0.61 

0.37 

Maximum  Deviation 
(deg) 

0.15 

0.25 

0.12 

0.63 

0.96 

0.58 

1.81 

3.38 

2.13 

worst-case  position  uncertainty  (0.6“)  would  cause  less  than  a one- 
half  degree  error  in  any  attitude  angle  estimate.  Case  II  included 
magnetometer  position  uncertainty  and  noise.  The  uncertainty  num- 
bers were  the  same  as  those  used  in  Case  I.  The  noise  was  Gaussian 
with  a one-sigma  rating  of  1 mA/m.  The  maximum  increase  in  the 
average  estimate  errors,  due  to  the  addition  of  the  noise,  was  only 
0.14*.  Case  III  was  the  same  as  Case  II,  except  that  the  input 
noise  had  a one-sigma  rating  of  5 mA/m.  Of  course  this  case  pro- 
duced the  largest  attitude  angle  estimate  errors.  However,  on  the 
average  the  roll,  pitch,  and  yaw  angle  estimates  remained  less  than 
1*  away  from  their  true  values. 
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6.  APPLICATION  OF  LEAST  SQUARES  TECHNIQUE 


Attitude  estimation  results  from  several  APL  spacecraft 
(Triad,  TIP-II,  and  TIP-lII)  are  presented  in  this  section.  Each 
satellite  was  designed  to  be  three-axis  attitude-stabilized.  Each 
was  configured  so  that  gravitational  and  gyroscopic  forces  would 
tend  to  align  the  satellite's  axes  with  the  local  vertical  system 
of  axes  (see  Fig.  1). 

The  estimated  roll,  pitch,  and  yaw  angles  express  quantita- 
tively the  angular  deviations  of  the  satellite's  axes  from  the  local 
vertical  set.  The  smaller  the  magnitudes  of  these  angles,  the  bet- 
ter the  attitude  stabilization. 


RESULTS  FROM  TRIAD 

The  Triad  satellite  (see  Fig.  5)  was  launched  into  a polar 
orbit  during  the  fall  of  1972  (Ref.  18).  Table  2 is  a chronology 
of  the  significant  events  prior  to  and  including  its  achievement 
of  three-axis  attitude  stabilization.  Included  in  the  table  are 
comments  that  point  out  the  anticipated  characteristics  of  the 
attitude  dynamics  that  are  compatible  with  each  event. 

The  attitude  angles  that  were  computed  according  to  the  algo- 
rithm outlined  earlier  are  listed  in  Tables  3 through  11.  Each 
table  contains  the  results  from  one  day's  collection  of  attitude 
data  during  passes  over  APL.  The  horizontal  lines  in  each  table 
separate  the  various  passes.  When  the  results  are  grouped  accord- 
ing to  the  events  listed  in  Table  2,  it  can  be  seen  that  the  atti- 
tude performance  was  as  expected.  The  best  attitude  stabilization 
was  observed  on  the  last  two  days,  271  and  272,  when  the  attitude 
dynamics  appeared  to  be  near  steady-state.  The  angle  between  the 
Z-axls  and  the  local  vertical  was  consistently  below  5®. 

For  each  attitude  estimation  result,  there  are  two  error  in- 
dicators. One  indicator  is  the  difference  between  the  theoretical 
and  observed  angles  between  the  geomegnetlc  field  and  sunline 


Ref.  18.  Space  Dept.  Staff  of  APL  and  the  Guidance  and  Con- 
trol Staff  of  Stanford  University,  "A  Satellite  Freed  of  All  but 
Gravitational  Forces:  TRIAD  I,"  Paper  No.  74-215,  presented  at 
AIAA  12th  Aerospace  Sciences  Mtg.,  Washington,  DC,  January  1974. 
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Table  2 

Chronology  of  Significant  Events  Affecting 
Triad  Attitude  Dynamics 


I  2 September  (Day  2A6) : Triad  launch 

Satellite  now  In  a spinning  and  tumbling  mode.  Roll, 
pitch,  and  yaw  angles  meaningless  at  this  time. 

II  4 September  (Day  248,  02  h 50  min):  magnetic  stabilization 
Initiated 

Z-coll  turned  on.  Momentum  wheel  revved  up  to  achieve 
0.25  slug-ft^  rpm  momentum.  Quality  of  stabilization 
can  be  determined  from  angles  In  third  column  from  left 
on  attitude  determination  printout.  Angle  of  180®  Is 
perfect  stabilization.  Momentum  wheel  should  provide 
yaw  control.  Roll,  pitch,  and  yaw  angles  are  now  mean- 
ingful during  passes  over  API. 

Ill  7 September  (Day  251,  15  h 13  min):  gravity-gradient  capture 

Satellite  booms  deployed.  Rotor  momentum  increased  to 
4.0  slug-ft2  rpm.  Amplitudes  of  roll,  pitch,  and  yaw 
angles  should  consistently  diminish  from  day  to  day. 

IV  14  September  (Day  258,  22  h) : trouble;  low-voltage  sensor 
switch  tripped 

Momentum  wheel  power  turned  off.  Its  angular  rate 
approaches  zero.  Yaw  axis  stabilization  Is  lost. 

V  15  September  (Day  259,  02  h 23  min):  recovery  period 

Momentum  wheel  power  Is  restored.  Large  pitch  angle 
libratlons  result  from  rotor  torque.  Large  roll  and 
yaw  libratlons  result  from  gyro  coupling  of  rotor 
momentum.  Roll,  pitch,  and  yaw  angle  amplitudes  should 
again  diminish  from  day  to  day. 

VI  29  September  (Day  273):  trouble;  8-bit  A/D  converter  lost 

Cannot  obtain  DSAD  and  magnetometer  data  In  the  usual 
manner.  Other  possible  means  cannot  be  readily  Imple- 
mented. Generation  of  roll,  pitch,  and  yaw  angle  esti- 
mates Is  thus  terminated. 
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Triad  Attitude  Determination  Results,  10  September  1972 
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vectors.  The  other  Is  the  difference  between  the  theoretical 
and  observed  geomagnetic  field  magnitudes.  In  general,  the 
smaller  the  differences,  the  more  accurate  the  attitude  angle 
estimates.  The  observed  and  theoretical  data  are  listed  In 
columns  7 through  10  (from  left  to  right)  In  each  table. 


RESULTS  FROM  TIP-II 


The  second  satellite,  TIP-Il  (see  Fig.  6),  was  launched  In 
1975.  It  failed  to  achieve  the  desired  structural  configuration 
(for  gravity-gradient  stabilization)  because  of  a malfunction  of 
Its  extendible  boom.  The  resulting  moment-of-lnertla  distribution 
was  unfavorable  In  that  attitude  stabilization  about  the  desired 
spacecraft  X-,  Y-,  and  Z-axes  was  not  possible. 

The  following  attitude  estimation  results  on  23  September 
1976  were  typical  of  the  satellite's  dynamics  for  more  than  2 
months.  On  this  day,  at  3 p.m.  EST,  approximately  30  min  of  atti- 
tude sensor  data  were  processed.  A summary  of  the  attitude  esti- 
mation results  Includes  the  following. 

1.  A total  of  259  sets  of  digitized  attitude  sensor  data 
was  received  before  the  satellite  entered  the  night  portion  of  the 
orbit.  More  than  80%  of  the  sets  passed  all  validity  checks. 
Figure  7 Is  a plot  of  the  sun  sensor  data.  The  i(;  angle  (ordi- 
nate axis)  Is  the  angle  between  the  sunllne  vector  and  the  vehicle 
Z-axls.  The  azimuth  angle  (abscissa  axis)  Is  the  angle  between 
the  X-Y  component  of  the  sunllne  vector  and  the  vehicle's  X-axls. 
The  time  history  of  the  magnetometer  data  Is  shown  In  Fig.  8. 

2.  The  excellent  agreement  between  the  computed  (theoreti- 
cal) and  measured  (observed)  magnetic  field  and  sun  data  Is  shown 
In  Figs.  9 and  10. 

3.  The  satellite  attitude  dynamics  was  determined  to  be  a 
combination  of  a 3-  to  5-revolutlon  per  orbit  pitch-axis  tumble 
and  a relatively  rapid  coning  of  the  pitch  axis  In  Inertial  space. 
Gravity-gradient  stabilization  of  the  reference  Z-axls,  as  Indi- 
cated by  the  pitch  tumble,  had  not  occurred.  Indeed  the  dynamics 
appeared  more  representative  of  a spin-stabilized  satellite.  Fig- 
ures 11  to  15  are  plots  of  the  estimated  attitude  angles.  The 
last  plot  shows  clearly  the  coning  motion  of  the  spacecraft's 
Y-axls  in  Inertial  space.  The  declination  and  right  ascension 
angles  are  referenced  to  the  geocentric  reference  system  of  axes 
(Z  Is  the  North  Pole,  X Is  the  first  point  of  Arles,  and  Y is  the 
vector-completing  right-hand  set). 
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Fig.  7 Plot  of  TIP-II  Sun  Sensor  Data  (23  September  1976) 
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Fig.  12  Plot  of  TIP-II  Pitch  Attitude  Angle  (23  September  1976) 
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Fig.  13  Plot  of  TIP-II  Yaw  Attitude  Angle  (23  September  1976) 
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Fig.  14  Plot  of  TIP-II  Roll,  Pitch,  and  Yaw  Attitude  Angles 
(23  September  1976) 
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RESULTS  FROM  TIP-III 


TIP-III  (the  same  configuration  as  TIP-II)  was  more  success- 
ful than  Its  twin  in  achieving  gravity-gradient  stabilization.  Dur 
Ing  the  early  part  of  March  1977,  the  satellite’s  boom  was  success- 
fully extended,  and  a right-slde-up  stabilization  was  achieved. 

The  results  presented  here  are  from  10  March  1977  (12  noon) , when 
the  attitude  was  monitored  over  a period  of  three  successive  orbits 
The  results  Included  the  following. 

1.  The  sun  sensor  data  are  displayed  as  time  history  plots 
in  Figs.  16  and  17.  The  i|)  and  azimuth  angles  are  as  defined  pre- 
viously. The  noon  (ip  is  close  to  zero)  and  midnight  (i|/  is  approxi- 
mately 180®)  portions  of  each  orbit  are  readily  discernible  in  the 
figures . 


2.  The  magnetometer  data  are  shown  in  Fig.  18.  The  Z-mag- 
netometer  shows  the  approximate  times  the  spacecraft  crossed  the 
North  and  South  Poles  with  its  minimum  (-450  mA/m)  and  maximum 
(+450  mA/m)  readings,  respectively. 

3.  The  peak  roll  angle  was  30*  (Fig.  19).  The  frequency  of 
roll  oscillation  was  approximately  2 revolutions  per  orbit. 

4.  The  peak  pitch  angle  was  15®  (Fig.  20).  The  frequency 
of  pitch  libratlon  was  approximately  1.7  revolution  per  orbit. 

5.  During  this  time,  the  constant  speed  rotor  was  off. 

The  satellite's  inertia  ellipsoid  was  very  much  like  that  of  a 
dumbbell.  Without  the  rotor,  the  moment-of-lnertla  distribution 
would  not  tend  to  stabilize  very  well  (if  at  all)  with  the  Y-axis 
aligned  with  the  orbit  normal.  Thus  it  was  not  surprising  when 
the  spacecraft,  on  this  day,  was  observed  performing  360®  rota- 
tions in  yaw  approximately  once  every  3 h (Fig.  21) . The  rotation 
was  negative. 
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Fig.  20  TIP-Ill  Pitch  Attitude  Angle  (10  March  1977) 
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7.  CONCLUSION 


A rather  complex  formulation  for  estimating  the  attitude 
of  a satellite  In  a least-squares  sense  has  been  Implemented  at 
APL.  The  technique  has  been  successfully  used  to  monitor  and 
assess  the  In-orblt  attitude  dynamics  of  several  gravlty-gradlent- 
stablllzed  APL  spacecraft.  Its  realization  represents  a signifi- 
cant addition  to  the  attitude  estimation  techniques  that  have  been 
developed  to  support  APL  satellite  operations  for  more  than  a decade. 
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Appendix  A 

DERIVATION  OF  LEAST-SQUARES  SOLUTION 

In  this  Appendix,  the  least-squares  solution  to  the  atti- 
tude estimation  problem  will  be  derived.  The  discussion  will 
parallel  the  development  in  Ref.  7.  Capital  alphabets  will  denote 
matrices,  and  small  alphabets  will  denote  column  vectors.  Scalars 
will  be  appropriately  identified. 

Given  two  sets  of  n unit  vectors  m^ , m_ , ...  , m and  m*, 
* * i z n 1 

m2,  ...  , m^,  where  n i 2,  find  the  rotation  matrix  A (l.e.  , the 

orthogonal  matrix  with  determinant  -bl)  that  brings  the  first  set 
into  the  best  least-squares  coincidence  with  the  second.  That  is, 
find  A,  which  minimizes 


n 

E "*■'^3 


The  problem  has  arisen  in  the  estimation  of  the  attitude  of 
a satellite  by  using  the  unit  vectors  (m*)  of  objects  as  observed 
in  a satellite  fixed  reference  system  and  the  unit  vectors  (m)  of 
the  same  objects  in  a known  reference  system.  A is  then  a least- 
squares  estimate  of  the  rotation  matrix  that  carries  the  known 
reference  system  into  the  satellite  fixed  reference  system. 

Let  k denote  the  number  of  elements  of  the  column  vectors 
m, , ...  , m , mf,  ...  , m and  let  M and  M*  denote  the  two  k 
1 " A A 

by  n matrices  whose  columns  are  m, , ...  , m and  m m , 

respectively.  ^ n 1 n’ 

Now  define  Q(A)  as  the  sum  of  squares  to  be  minimized,  l.e.. 
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where  tr  denotes  the  trace  function  and  a superscript  T denotes 
transposition. 


Q(A)  can  be  expanded  and  rewritten  as 


Q(A) 


tr 


T T * 

M A ) (M  - AM) 

2tr(M^AV). 


* 


tr(M  M ) + tr(M  M) 


Since  the  first  two  terms  are  Independent  of  A,  Q(A)  Is  minimized 

T T * 

by  maximizing  F(A)  - tr(M  AM).  This  function,  F(A) , may  be 
written  as  (due  to  the  preservation  of  commutativity  under  the 
trace  operation) 


F(A)  - tr(A^  M*  M^)  . 


It  Is  well  known  that  an  arbitrary  real-square  matrix,  P, 
can  be  written  as  a matrix  product,  US,  where  U Is  orthogonal 
and  S Is  symmetric  and  positive  semldeflnlte . Furthermore,  If 
P Is  nonsingular,  U Is  uniquely  defined  and  S Is  positive  defi- 
nite. If  P Is  singular,  U Is  not  unique  and  the  problem  Is 
essentially  Indeterminate. 

* T T 

Applying  the  above  to  P ■ M M , we  have  F(A)  = tr(A  US). 

Since  S Is  symmetric,  there  Is  an  orthogonal  matrix  G such 

that  (CSC'^)  Is  a diagonal  matrix  D,  whose  diagonal  elements  d]^, 

...  , d.  are  arranged  In  decreasing  order.  All  d,  are  non-nega- 

^ J T T 

tlve  since  S Is  positive  semldeflnlte.  Now  letting  X “ (GA  UG  ) , 

we  obtain 

k 

F(A)  - trCA^G^DG)  - tr(GA^UG^D)  - tr(XD)  - d^x^^  . 

1-1 

Since  F(A)  la  a linear  function  of  the  non-negative  numbers 
d^,  ...  , d^,  its  maximum  Is  attained  when  the  diagonal  elements 

of  X attain  their  maximum  values.  Because  X is  an  orthogonal 
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matrix,  all  of  Its  elements  are  between  -1  and  +1  In  value.  Thus 
F(A)  Is  maximized  when  ■ +1  and  x^j  = 0 for  1 j . 

T T 

Because  det  A is  required  to  be  +1,  the  det  X “ det  (GA  UG  ) 
- (det  G)2(det  A)  (det  U)  - det  U.  If  det  U »»  -1,  it  is  required 
that  det  X = -1  and  it  is  not  difficult  to  see  that 


= I - 2H 


is  a solution  (since  is  d2  i ^ matrix,  I,  is  the 

kth-order  identity  matrix  and  ^ is  the  k-lst  order  identity 

matrix.  H is  a kth-order  square  matrix  with  +1  as  the  (k,  k)  ele- 
ment and  zero  everywhere  else.  For  the  case  where  det  U = +1, 

X = I.  Now  letting  Xq  be  the  matrix  that  maximizes  F(A)  (either 
X=IorX*I-2H,  according  to  det  U = +1  or  -1),  then 


T T T T 

Xq  = G Aq  U G , or  Aq  = U G Xq  G 


is  a rotation  matrix  that  minimizes  the  sum  of  squares  Q(A).  If 

* T 

P - M M is  nonsingular,  it  is  the  unique  rotation  matrix  that 
does  so. 

In  the  solution  developed  by  R.  H.  Wessner  (Hughes  Aircraft 

* T 

Co.),  he  points  out  that  if  det  P/0,  P»MM  “US,  and 


T -1  T *5  T *5 

U - (P  ) (P  P)  , S - (P  P)  . 

T <j' 

where  (P  P)  is  the  symmetric  square  root  of  (P  P)  with  positive 
eigenvalues . 

Hence  the  solutions  for  A are 


^ * T 


Aq  - U - (MM  ) (MM  MM)  for  det  P > 0,  and 
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T -1 


(MM  ) 


T *J 

at 

(MM  M M ) (I 


2C  HG)  for  det  P < 0. 


The  matrix  G 
P^P. 


Is  the  modal  matrix  of  eigenvectors  for  the  matrix 
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Appendix  B 

ROTATIONAL  TRANSFORMATION  BETWEEN  LOCAL 
VERTICAL  AND  GEOCENTRIC  REFERENCE  SYSTEMS 


The  mathematical  formulations  used  In  computing  the  sunllne 
vector  and  geomagnetic  field  vector  are  referenced  to  a geocentric 
coordinate  system  of  axes  (Z  Is  the  North  Pole,  X Is  the  first 
point  of  Arles,  and  Y Is  the  vector  that  completes  the  right-hand 
set).  This  means  that  the  vectors  are  In  geocentric  coordinates 
rather  than  the  desired  local  vertical  coordinates.  One  way  to 
solve  this  problem  Is  to  compute  an  orthogonal  transformation  that 
defines  the  orientation  between  the  two  coordinate  systems. 

The  derivation  of  this  matrix,  denoted  £,  begins  with  the 
observation  that  the  local  vertical  system  of  axes  (X^,  Yj^,  and 
Zf^)  are  parallel  with  orbit  axes  whose  direction  cosines  are  a 
function  of  the  standard  Kepler  orbit  elements,  fl,  u),  1,  and  f (de- 
fined below).  These  orbit  axes,  which  will  be  called  Yq,  and 
are  Identified  next. 

Figure  B-1  shows  the  geometry  of  the  satellite's  orbit,  as 
well  as  the  location  of  the  3^,  Y^,  and  ^ axes.  The  sequence  of 

angular  rotations  from  the  geocentric  coordinate  axes  to  the  set 
of  orbit  axes  Is  as  follows: 

1.  A rotation,  , about  the  axis  Is  performed.  The 

angle,  S2,  Is  usually  called  the  longitude  of  the  ascending  mode. 

2.  A rotation,  ^(1),  about  the  new  X-axls  (also  called 

the  line  of  nodes)  Is  next  performed.  The  angle,  1,  Is  the  orbit 
Inclination. 

3.  A rotation,  , about  the  new  Z-axls  is  next  per- 

formed. The  angle,  u,  Is  called  the  argument  of  perigee.  The  new 
X-axls  (line  of  apsides)  of  this  system  defines  the  points  In  orbit 
of  the  satellite's  closest  (perigee)  and  farthest  (apogee)  ap- 
proaches to  the  earth's  mass  center. 

4.  The  last  rotation,  ^^(f),  about  the  new  Z-axls  Is  per- 
formed. The  angle,  f.  Is  called  the  true  anomaly.  The  X-axls  of 
this  system  Intersects  the  orbit  path  at  the  point  where  the  satel- 
lite's center  of  mass  Is  situated. 
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North  Pole 


Orbit  normal 


Line  of 
apsides 


(first  p>oint 
of  Aries) 


Earth  equa^ 


Line  of  nodes 


Orientation  of  Local  Vertical  System  Relative  to  Geocentric 
Reference  System 
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The  product  of  these  rotations  is  a direction  cosine  matrix, 
D,  that  can  be  written  as 

m ’ 


D - R^(f)R2(a))Rj(i)R^(n) 
- R3(F)R^(i)R^(P)  , 


where  the  last  two  Z-axls  rotations  are  equivalent  to  a single 
Z-axis  rotation  of  magnitude,  F ■ f + oi.  Note  that  the  axis,  3Cq, 
is  the  outbound  local  vertical  and  that  2^  is  the  orbit  normal. 

Identification  of  the  local  vertical  axes,  Xj,,  Yj^,  and  ^ is  thus 
equivalent  to  a redefinition  of  the  2^,  and  ^ axes.  This  can 

be  accomplished  using  two  90*  rotations.  The  first,  22(90*), 
places  a Z-axis  along  the  direction  of  3^.  The  second  rotation, 
23(90*),  about  the  new  Z-axis  places  the  new  Y-axis  along  the  di- 
rection of  Thus  the  total  transformation  matrix,  C,  from  the 

geocentric  reference  system  to  the  local  vertical  reference  sys- 
tem is 


C - R,(90*)R.,(90*)D 

- -3  “2  ■ 


or 


C 


-fJ  F 

- fi  F 1 

fi  1 

fi  F 

- n F 1 V 

c s 

see 

s s 

c c 

s s c \ 

-n  F 

+ r.  F 1 

-n  1 

n F 

+ fi  F 1 1 . 

s a 

c c c 

c s 

s c 

c s c 1 

F 1 

1 

F 1 

/ 

c s 

c 

s s 

/ 

The  determination  of  the  orbit  elements  Is  an  orbit  deter- 
mination problem.  In  the  attitude  estimation  scheme,  these  are 
assumed  known. 
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NOMENCLATURE 


- Orthogonal  transformation  matrix  from  a local  vertical 
reference  system  to  the  reference  axes  fixed  in  the 
satellite 

■ Orthogonal  transformation  matrix  from  the  geocentric 
coordinate  system  to  the  local  vertical  reference  system 

“ Orbit  longitude  of  the  ascending  node 

« Orbit  argument  of  perigee 

“ Orbit  inclination 

“ Orbit  true  anomaly 


pedal  Notation 


1 


(a) 


0 


a 

c 


-a 

s 


0 


••  matrix  representation  of  a positive 
rotation  of  a radians  about  a 1- 
or  X-axls 


R2(a) 


0 

1 

0 


* matrix  representation  of  a positive 
rotation  of  a radians  about  a 2- 
or  Y-axls 


a ■ Sine  of  a 

s 


a 

c 


0 

0 

1 


■ matrix  representation  of  a positive 
rotation  of  a radians  about  a 3- 
or  Z-axls 


Cosine  of  a 


1C  le 


tm  JOUHS  HOPHINS  UN»V€«S>Ty 

APPLIED  PHYSICS  LABORATORY 

l«4in(k  MAHviANf) 


U 


1 


- Column  vector,  u,  In  the  satellite  coordinate  system 

■ Column  vector,  In  the  geocentric  coordinate  system 

■ Column  vector,  u.  In  the  local  vertical  coordinate  system 

- Column  vector,  u,  in  the  orbit  coordinate  system 

■ Transpose  of  the  matrix,  U (applies  to  column  vectors 
also) 

■ Inverse  of  the  square  matrix,  U 


k 


-z 

J-1 


a 


tr  U 


trace  function  of  the  kth  order  matrix,  U 
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